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1 . 0 INTRODUCTION 


This project is the initial stage of a study to determine improved values 
and definitions to be used for thermal environmental design parameters for 
a spacecraft in near-earth orbit. The current criteria is inadequate due 
to incorrectly assumed confidence limits, the lack of interrelation between 
albedo and earth thermal radiation, and no altitude dependence. It is 
planned at a future time to develop the improved criteria through a care- 
ful statistical study of known data. 

As a starting point, this project was established to convert existing 
satellite data from selected orbits for which both instantaneous values 
of the earth thermal and albedo radiation could be obtained, into data 
which could be easily handled for the planned computer oriented statistic 
studies . 

The Medium Resolution Infrared Radiometer (MRIR) carried onboard the NIMBUS II 
spacecraft seemed particularly suited to this study since it contained both 
a 5-30 micron channel to provide a measurement of the earth thermal radiation 
and a .2-4 micron channel to provide the radiance of the solar reflected 
radiation (albedo). However, because the lifetime of this instrument was 
so short it was necessary to use as much of its limited information as poss- 
ible, to develop a method of utilizing the data from the later NIMBUS III 
MRIR. The NIMBUS III MRIR sensor had five channels for observing earth 
radiation. 
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Longwave 
Albedo 

no single channel covers the entire 
II MRIR sensor. 

Having conducted a thorough literature search in this field, it was found 
that the NASA Technical Note #D-7249 "The Radiation Balance of the Earth- 
Atmosphere System from Nimbus 3 Radiation Measurements" by Raschke, 

Vonder Haar, Pasternak, and Bandeen was most applicable to this project. 
Goddard Space Flight Center supplied a listing of the program used by 
Raschke, Vonder Haar et al. to generate their data and the techniques con- 
tained within this program were used as a basis for developing the programs 
used in this project. 

A basic concept of the GSFC study was to develop an algorithm whose co- 
efficients could be used to derive a total earth thermal radiation based 
on a mathematical relationship obtained from studying 160 different 
atmospheric models. This replaces the need to develop a single algorithm 
by performing a strictly mathematical correlation of the various spectral 
bands to the broad band channel of the NIMBUS II data. To produce an out- 
put compatible with other studies performed in this field, this same 
technique was used for this project. 



0.2 


4.0 


It can be seen, unfortunately, that 
longwave spectrum as did the Nimbus 
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The end product produced by this project is 100 albedo and 100 earth 
thermal radiation grid maps on 7 track digital magnetic tape. A grid 
map is a matrix of elements where location in the matrix is related 
to an earth latitude and longitude. Each value in the matrix is the 
average of the data over the matrix element area. The size of the 
matrix element area varies depending on the latitude being observed. 

The area is 1° latitude by 1° longitude at latitudes less than 60°; 

1° latitude by 2° longitude at latitudes between 60° to 80°, and one 
area for each polar region at latitudes of 80° and over. Each map con- 
tains the values obtained during a 24 hour period over the entire earth. 

No extrapolation was performed on the data, since this will allow greater 
flexibility for later simulation and Time Average Effective Value (TAEV) 
computations. The confidence interval widths provided with the statis- 
tics for each matrix element area of the earth's surface (data bin) can 
also be used at some future time for dynamic smoothing of the bin data. 

It is thought that such autoregressive techniques are superior to any 
interpolation/extrapolation schemes which could permanently alter the 
values contained in the constructed data base. 

2.0 BRIEF SURVEY OF RESULTS 

A brief sample of the output statistics is plotted in the following graphs 
(Figures 2.0-1 and 2.0-2) to illustrate the type of data generated by the 
program. The information plotted is the bin averages. 
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Figure 2.0-1 

Outgoing Longwave Radiation (cal/cm min) 
Measured by NIMBUS III Day 127 , Orbits 315-317 
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Figure 2.0-2 

Albedo (Percent) Daily Averages 
Measured by NIMBUS III Day 127, Orbits 315-317 
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Figure 2.0-2 (cont.) 

Albedo (percent) Daily Averages 
Measured by NIMBUS III Day 127, Orbits 315-317 
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A check on the validity of the data produced by the program was made 
by comparing the data with contour plots of albedo and longwave radia- 
tion contained in Appendix D of NASA Technical Note NASA TN D-7249. 

The values from the contour plots used for the comparison are approxi- 
mate estimates, however, they do provide indication that the computed 
values are comparable . 

The information is summarized in the following two tables: 


Lat . 



Long. Range 

Computed 

Range 

Data Obtained 
From 

NASA TN D-7249 

-45 

90-130 

. 19-. 32 

.30-. 33 

-40 

90-130 

.17-. 32 

.33-. 36 

j -30 

90-130 

.18-. 37 

.33-. 36 

-20/-18 

83-118 

.37-. 42 

.39-. 42 

-10 

100-128 

.40-. 44 

.39-. 42 

1 0 

102-150 

. 23-. 43 

.33-. 36 

10 

! 106-142 

.14-. 40 

. 36-. 39 


Table 2-1 

Outgoing Longwave. Nimbus III. Orbits 315-317 





Data Obtained 



Computed 

From 

1 IPV1 

Long, Range 

Range % 

NASA TD D-7249 

-58 

80-126 

44-59 

50 

-44 

28-134 

33-60 

40 

-20 

96-134 

5-36 

20 

1 

104-142 

9-19 

20 

10 

106-142 

6-68 

20-30 


Table 2-2 

Albedo (Percent), Nimbus III. Orbits 315-317 
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3.0 data processing program 


The program developed for this project, computes estimates of albedo and 
longwave radiation by extracting the MRIR data from NIMBUS Meteorological 
Radiation Tapes. 

This data is stored by "swaths", each representing a sweep of the rotat- 
ing scanner. The location of the satellite, time and scan angle associated 
with the data are sufficient to enable the location of the earth reference 
area observed by the sensors to be determined. 

Following the processing of the data the summary statistics of the radia- 
tion estimates are accumulated onto the digital magnetic output tape. 

These statistics are detailed in Section 5.2 . 

A major functional flow of the program is presented in Section 3.1 
followed by a detailed description of the steps involved in generating the 
estimates . 
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3.1 FUNCTIONAL FLOW DIAGRAM 


° Set file counter 
to one 


° Numoer of swaths/ 
record 

° Size parameters 
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Program SWATH 


LOOP FOR 
NUMBER OF 
DATA POINTS 



DETERMINE 
LIMITS FOR 
DATA POINTS 
f (NADIR) 


* 



DETERMINE 
BIN NUMBER 
LAT/LONG 
INTERPOLATION 





CALL BALANCE 
TO COMPUTE 
PRIMARY OUTPUT 
QUANTITIES 

- ' 

1 _ 


DETERMINE 
NIGHT/ DAY 
CONDITION & 
ACCUMULATE 
BIN 

STATISTICS 



i 


II 





3.1.1 Calculation of the Total Outgoing Longwave Radiation 

The two primary processing steps performed by the program in computing, the 
total outgoing longwave radiation are as follows: 

Step (a) Estimate the tota l Longwave Radiation for the entire bandwidth^ 
from a NIMBUS II derived correlation formula, 

The measured data consists of the filtered radiance observed concurrently 
by the four channels of the MRIR instrument. Assume that the values 
from each channel are N^, N 2 , N 3 and N 4 respectively* 


KL "■ ®(c '' 'Xi Nj 


M * + 


where 

<X 0 = -.35651 X10~ 2 

= .37442 XlCT 2 

= .30818 X10' 3 

0^ = -.92471 X10" 5 

= .013466 

0(i = .20913 X10- 2 

<*£ = .3051 X10“ 2 


Step (b) Compute the Outgoing Flux Density (N). 

This is obtained by applying a correction which is dependent on satellite 
location (reference Appendix A for a definition of terms)* 
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The outgoing flux is based on the following estimator: 

~ ZTT / M t S\^e J B ~ 1 * V 

r> 

'•'/t 


where 


where 


^ / £('8j fosB riKN ^ 
0 

^ S') 


f re") = 


ru-co) 


Uo,i:to 






^ 5 ~'io ( 0 - o ^ 


z 3 

- \ r-fi.B-t fB z ft I- fis® 


This polynomial model is derived from NIMBUS II data and represents the 
drop in emitted radiation as 9 increases. This radiation model is 
referred to as " 1 imb- darken ing" . 


It should be noted that the final outgoing longwave flux represents an. 
estimate of the total radiation. 


The normal radiation ( N 5 _ 30 (#=0)) is also derived from NIMBUS II 
measurements . 


Computationally the program uses the following: 

Forj^^ 60° and; _ 

© <10°, N = 3.05 N t 
10 Q <e<20 o , N = 3.054 N t/ 985 

20°<a<30°, = 3.054 N t/>975 


For 60 u and; 


05= 30°, FT = 3.17 N t 
0> 30°, N = 3.17 N t y >995 
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3.1.2 Calculation of the Incoming/Out going Solar Radiation (Albedo) 


For shortwave or albedo radiation the program computes the average of 
the outgoing flux densities from the observed data. 

The measured data (N ) consists of the data from the 0. 2-4.0yuj band 
converted to units of watts/ (master .) . The primary processing steps 
performed by the program in computing the albedo radiation are as follows 
Step (a) C ompu t e the Bidirection a l Reflectance . 


The Bidirectional Reflectance is the ratio of the measured data emerging 
from a surface element in the direction of the satellite to the total 


radiation entering normal to the surface from space. 


satellite 

location 




•InT 


zenith 
angle of 
measurement 


4 A 


azimuth angle 
of measurement 


solar zenith 
angle 


4 4 


earth reference 
longitude 


earth reference 
latitude 


<> a L coit' 


time of 

measurement 


N - N( ^ ^ j <$> t ) watts/m^ster 

measured data from the,2-4^ band 


So 

L 


- 860.58 watts/m^ filtered extra terrestrial 
radiance of sun (0,2 - 4.0//) 



mean sun-to-earth distance 
true sun-to-earth distance 


★"Primes" denotes measured data. 
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Step (b) C ompu te the appr oximat e directional reflectance, 

The directional reflectance can be obtained by using the following 
theoretical equation: 

2 tt yt 

r(^) - j j %) si'ab cos 9 ciSd p 

e a 


However, instead of integration, each observed bidirectional reflect- 
ance data point can be converted to a directional reflectance point 
at a specific surface element by using empirical or measured data. The 

goal being to remove directional (F,^ ) dependence by use of the ratio 

_ 

which is available as published tabular data. 

The data is stored by: 

(1) Surface (cloud-land, ocean, etc.) 

( 2 ) % 

(3) B 

(4) p 

Definition of item 1 is obtained from the range of p / 

e.g. 

(1) Ocean TT P < .10 and Tfe, 

(Tb oeing the concurrently observed blackbody temperature) 

(2) Snow | > <b5 6 and 

(3) Cloud/land (anything else)t 

/ 

Directional reflectance is computed by 


rm 

r_Ti.T ”1 

jtUc) 

[_ ir y) 

Table F- 1 

TftWe F-^ 
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It can be seen that r is not a single value, but is a function of two 
tables; one a function of the surviving independent variable £ and the 

/ 

other a function of the observed or concomitant angular values >. 

The tables are contained in Appendix F. 

Step (c) Est imate the Inco min g/ Out going Radiation 

The estimates are derived utilizing the following 

Incoming (t) = (S 0 L) f ol 1= ("t ) 

Outgoing (t) = rostCt') 

\ = solar constant 

! L = reciprocal of sun-earth distance squared in astronomical units 

These are integrated over the daylight period for day 1 t ' to arrive at an 
estimate of the daily averages (see Appendix E). 
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3.2 THE NIMBUS III DAYS PROCESSED 

The following table contains the day numbers for the NIMBUS III data that 
was processed. An attempt was made to utilize equal daily increments, 
however this was tempered by determining which days had more representa- 
tive data and the availability of data tapes. 


1969 : 1970 r 


DAY NUMBER 

DAY NO. 

! 

} 107 

163 

222 

21 

I 109 

164 

223 

23 

: no 

166 

224 

25 

112 

168 

226 

28 

114 

170 

228 

29 

116 

172 

230 

32 

118 

174 

232 

34 

120 

176 

234 


122 

178 

236 


124 

180 

238 


126 

182 

239 


128 

184 

241 1 


130 

186 

244 i 


132 

188 

246 S 


134 

190 

249 


136 

192 

252 j 


138 

194 

256 


140 

196 

258 


142 

198 

261 


144 

201 

264 


147 

202 

267 


148 

204 

270 


149 

206 

273 


151 

208 

277 


152 

210 

281 


154 

212 

285 


156 

214 

289 


158 

216 

293 


160 

218 

298 


161 

i I 

220 I 
i 

303 



Table 3-1 Day Numbers of Nimbus III Data Processed 
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4.0 INPUT DATA STRUCTURE 


The input data form and content is described in detail in the "NIMBUS III 
User's Guide 11 prepared by the NIMBUS Project, Goddard Space Flight Center, 
National Aeronautics and Space Administration. Section 4 (PP 67-107) 
deals with the MRIR data which is the principal source of data. The for- 
mat of the NMRT-MRIR tape is described on PP 99-103 of that document. 

5.0 OUTPUT DATA FORMATS 

The program produces a digital magnetic tape output containing the daily 
averages for each bin and selectable line printer output of either 
summary or detail bin statistics. 

The digital magnetic tape is the normal output of the program, the printer 
output being used during the checkout of the program. However printouts 
can be requested if desired. 


5.1 DIGITAL MAGNETIC TAPE FORMAT 

The output tape is a standard 7 track, 800 bpi digital magnetic tape. 
Each tape may contain more than one data file with each file containing 
the compiled statistics for one satellite day. The magnetic tape is in 
Univac 1108 compatible binary format (e.e., formatted 36 bits/word). 


1 1 

1 


e" 

j E 

DATA FILE # 1 

:0DATA FILE #20 / DATA FILE #N 

0 

lo 



F 

iF 

Fig. 5.J 

7 / 

L-l Output Tape Gross Format 
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TYPE 1' 

TYPE 2 (TYPE 2“ 

fTl 

TYPE 2 

TYPE 2 i TYPE 3 

, 

E> 

RECORD 

RECORD j RECORD 


RECORD 

RECORD i RECORD 

0 


#1 | #2 

—LX 

#251 

#252 I 1 

' F 


Fig, 5,1-2 Data File Gross Format 


As Figure 5.1-2 depicts there are three different types of data records 
in each file. Every record is equal in size and consists of 1200 words 
(36 bits per word). 


Type 1 

The first record of each file is type 1. Word 1 contains 
the value of the satellite day processed and is the only 
information contained in this record. 


Type 2 

Records 2 through 253 of each file are type 2. Each type 2 
record contains the statistics for 200 bins, formatted as 
6 words for each bin. 

The format of each of the 6 word groups is as follows: 


woSr^Pos 
no. — ^ 

35 18 

17 . 0 

1 

Incoming Mean 

Outgoing Mean 

2 

Nighttime Mean 

Daytime Mean 

3 

Outgoing Standard 
Deviation 

Nighttime Standard 
Deviation 

4 

Daytime Standard 
Deviation 

Outgoing 95% 
Confidence 

5 

Nighttime 95% 
Confidence 

Daytime 95% 
Confidence 

6 

Incoming /Outgo ing 
Population 

■ Nighttime 
Population 

Daytime 

Population 


Pofl 

35 24 

23 12 

11 0 


Table 5-1 

Summary Statistics Format 

NOTE: All the data values contained in words 1 through 5 

are scaled by 218; and the value of incoming and 
outgoing populations are assumed to be equal. 
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Type 3 

The last or record 254 of each file is a type 3 record 
which contains statistics for the north and south polar 
regions. The first 6 words of this record relate to the 
arctic region; the second 6 words, to the antarctic. The 
13th word contains information for both polar regions. 

In each case words 1-5 are identical in form to words 
1-5 of the type 2 record. 

Words 6 and 12 of a type 3 record have the following configuration. 


BIT POSITION 

35 18 

17 0 

OUTGOING POPULATION 

NIGHTTIME POPULATION 


Word 13 of the type 3 record has the following configuration: 


BIT POSITION 


35 18 

r - — - 

17 0 1 

NORTH POLE DAYTIME 
POPULATION 

SOUTH POLE DAYTIME 
! POPULATION i 


Records 2-253 of each file each contain the statistics for 200 bins, for 
a total of 50400 bins. Record 2 contains bins numbered 1-200, Record 3 
contains those numbered 201-400, etc. 
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To relate bin number to the latitude, longitude of the bin, consider 
that the earth is divided into 5 areas: 


Area 1 - Anarctic - Latitude < -80°. Total 1 bin 

Area 2 - Latitude > -80° but < -60°. This region 

is divided into 1° latitude by 2° longitude, 
regions. Total 3600 bins. 

Area 3 - Latitude 5 -60° but < +60°. This region is 

divided into 1° latitude by 1° longitude regions. 
Total 43200 bins. 

Area 4 - Latitude ^>+60° but < 80°. This region is 

divided into 1° latitude by 2° longitude regions. 
Total 3600 bins. 

Area 5 - Arctic - Latitude > 80°. Total 1 bin. 


Records 2 through 19 contain the data gathered from Area 2. 

Records 20 through 235 contain the data gathered from Area 3. 

Records 236 through 253 contain the data gathered from Area 4. 

Record 254 contains the data gathered from Areas 1 and 5. 

For example, the statistics for the bins from Area 2 are formatted on the 
tape as shown in Table 5-2 


RECORD 

6 WORD 
BLOCK # 


ES9 


2 

1 

l 

-80° 

0 D -Long. 

< 2° 

2 

2 

2 

•80° 

Long. 

-=4° 

2 

200 

200 

-79° 

38°5 Long, 

«r 40° 

3 

1 

201 


40°f Long. 

42° 

3 

2 

202 

■as 

42°f Long. 

44° 

3 

200 

400 

-78° 

78°^ Long. 

80° 

19 

1 

3401 

-60° 

320°^ Long. 

322° 

19 

2 

3402 

“60° 

322°^ Long. 

«= 324° 

* 

• 

; 

' 



19 

200 

3600 

-61° 

358°« Ung. 

360° 


Table 5r2 Bin~Record Layout Example 
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5.2 LINE PRINTER OUTPUT 


The following items are contained within the statistical printout from 
the program. Figure 5.2-1 is a sample printout. 

1. Daily average of incoming solar radiation (cal/ (cm 2 min)) 

2. Daily average of reflected solar radiation (cal/ cm 2 min)) 

3. Nighttime outgoing longwave radiation (cal/ (cm 2 min)) 

4. Daytime outgoing longwave radiation (cal/(cm^ min)) 

5. Population count for 1-4 

6. Sample standard deviation for (2) 

7. Sample standard deviation for (3) 

8. Sample standard deviation for (4) 

9. 95% confidence band(^ width) for mean of (2) 

10. 95% confidence band (^ width) for mean of (3) 

11. 95% confidence band (^ width) for mean of (4) 

Some defining formulas are: 


1. Population count n 


2 . Mean 


-E«, 




3. Standard deviation S / (Xj - X)' 


n-1 


4. 100CiC% confidence bands for theoretical bin mean (yu ) 


i Width . 


where "t" is a (1 - y) percentile of the t distribution. 
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?9 
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•— -3ft 

0 . ■? m * 0 

o.orom 

0. 175R6 

31 
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3? 
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33 
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Appendix A 


SATELLITE -EARTH-SUN GEOMETRY 



P Sc^fU'-k Lo^tiotf 


Sun 



Greenwich hour angle of sun (measured westward) 
% Declination of sun 


Satellite 


f /\ 5 Longitude of subsatellite point 
tp s Latitude of subsatellite point 


Observed 

area 



Longitude of observed (reference) area 
Latitude of observed (reference) area 


R earth mean radius (6371km) 
H satellite altitude (1120km) 
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Appendix B 


LOCAL SUN -EARTH GEOMETRY 



"shaded 11 region 
represents the 
element -sun 
and 

element -zenith 
plane 


0 - observed area element 

(bin) centered at ("X,cp). 




9 


longitude of observed area 

latitude of observed area 

azimuth angle of measurement 

0 angle between two planes 

zenith-satellite and zenith-sun 
measured from the observed area. 

zenith angle of measurement 


a 


° angle between local zenith and satellite 
measured from the observed area. 

solar zenith angle 


0 angle between local zenith and sun measured 
from the observed area. 
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The actual observed data is N((9 j 0,^), which is the radiance emerg- 
ing from a surface element in the direction of a satellite located 
at a polar angle 6 and azimuth with respect to the sun-vertical 
plane (see figure). The sun is located at a polar angled with respect 
to the local vertical. 


The total radiation emanating from the surface element can be expressed 
as: W - f / 

From the diagram ip is measured from y. 


To convert to the familiar polar coordinates 

X = sinS cos ) - - sing sin^ 

y = sing sin (^-(-</>)) = sin& cos ip 
dxdy = - J ) d6> dc/J 


3* 

SB 

3% 


-cos^ sini p 

-sin 6 

cosfi 

3s> 

5St 

9P 


cos © cos^i 

-sing 

sinp 


sin# cos & 


Therefore, the total emerging radiation is 

ZTT Iffc 

«- // N (0j ^ ) sin # cos 0 d 0 dip 


O O 

The bidirectional reflectance, a unitless quantity less than one, is 
the ratio of N to the total energy entering normal to the surface from 
space 


N 

* ( > = SoL cost, 
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The directional reflectance is obtained by integrating out the 
directional dependence on spacecraft position 

r (£) = / / £) £u\P COS&d 0^ 0 

O O 

so that only the sun polar angle remains. 

If the surface is considered to be completely diffusing ^ ^ p ) 
then: 

nr u/ a 

\-C?) = / J p(%) «'*© ros9 d(9 d(J) - TT ^' ( ^ ) 

a o 
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Appendix C 


MATHEMATICAL RELATIONSHIPS 


C - 1 Solar Zenith Angle t, 

COS $ = X'S" = COSg C0$4> COSC^-^G.) +• SW S Sl'ocp 


^ ^ arc cos 8 cos <9 cosCV'Xo) + SmS J 

where /) and cp are the latitude and longitude of a surface element, /, 
is the Greenwich hour angle and $ is the declination of the 
sun. 


C - 2 Zenith Angle of Measurement 8 


SlYc^S Str. (H- 9^1 

r ' = rTh 

Sic\(TT-S) : SiynTT CdZ(-O) 
+• <osTr sm(-p) 

SIA0 








I-t & 7 - 50 cits not- pTOC^SS 



where is the nadir angle, \r the earth's mean radius, and h 
the satellite height 
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C - 3 Azimuth Angle of Measurement Ip 


<7 


( £> x X ) perpendicular to sun-zenith plane 
( X * ^0 perpendicular to zenith-satellite plane 

cos tp = cxx y >(s»y) _ 

| X * F*j j S x *| 

simplifying the numerator - 
(XX F ) . (SxX) = X.(F R xS x X) 

= X • T - ('Fr.^'x] 

= x • [ ) F*1 ( °3Q "s - l F*| cosl -1 x [j 
= iF R | [cos© co^- cost] 

where P is the angle between sun and satellite measured at ). 

| ”x x Fft 1 = \ Fr | snn© 

1 "S x ~X ] - S S 

Therefore 

COS e cosS - cos P 
sm ® sibi t, 


| x x Fr | \ S x X ! 


Lp : arrc.cS 


r cos d cos 
L S'VXft 


S — C O 
Sn\^ 



li'e-Cime. in — 

'Tis o-C subsate ll tie pcrmi: 

latitude of subset Unite poiwfc 
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P-_ 


Fr * s 
I P*\ 


S ^[KC^ccC^s-V)-^ cO-~XO!+^t^ SU ? 3 -^] 
u J V\-e'f<z C ar^& "s" c*te. o i*d e*jyrC*Si o v 'J 

v o v ’’ ros ! we" «** *> ' S rt>^ ' re Sy-e c V, veta . 


~ L k cc ^cC (>-->&; - 4-[k ce p s s(^ 5 -> g ) _ ^ SO ->^~\ 

+- [ k- S, 4 > a - S<*> 1 ^ 


3 k Ccp s C. ('■>s~^c 9 } +C 4 c O. ->o) •- gKC^'sCcp cC^^'j 2 (?)-}&) 

+ kV^sa^>6) hdVs z a-le»)-8kccp s ccp SOfeOte') s(M<b; 


^ k £ s'd) -eKS^s.S'? 


^ r 

= K CdP s + c z cp -£KC4>s Ccp[cC>S->&) ctW') + s()s-)c,)s(» fa )] 

+■ P ? s^ + - 2 K sc^scp 


~ K ^~ 1 “Z K^GCp ff CcpcC7s-’)la -f S^^sd? J 


z K 


? +) ~ 2 k [ccp.ccp c-O -^0 + Sc? s s 4 >] 


| =1 y/ K ? yl - 2 K feasts ro^q> Cas('}-‘}s') + 5 ia^c 





Appendix D 


SUN ANGLE CORRECTION 


The sun angle is corrected for refraction of solar radiation at various 
angles. This correction reduces the mathematically computed solar 
zenith angle. The table D-l contains the corrections used by the pro- 
gram: 

^ corrected = £ uncorrected - A 


i l. _ r T . - ' 

Range 

1 s 

Correction (/\) 

8B°5 

.35 

87 ?5 <■ t> < 88 ?5 

.3 

86?5<|>< 87?5 

.235 

85?5<^< 86°5 

.20 

! 83?5<S<.85?5 

.15 

J 80° 0 < 83 ?5 

.10 

? < 80?0 

0 


Table D-l Sun Angle Corrections 
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Appendix E 


NUMERICAL INTEGRATION TECHNIQUES 

The daily or 24 average of incoming solar radiation over a specific 
area Q,c$) is 

suv\ set 

WsCMjt'U-^ f S 0 L<-ttS^Ct*)<K*' 

where the quantities . L } ^ have been defined earlier* 

The flux density of solar radiation leaving the earth-atmosphere system 
over a specific area ) is 

= TZ J ros5ft*)at^ 

5'jY>Y\S£. 

where in both cases, t represents the day of the observation. 

Looking first at the 24 hour average of the incoming solar radiation. 

The theoretical expression is f-\ s > cp ) ~t > ) &S dlfSCVt bed OiboUti , 

where ’So = b'tb /{CXT\^ - fn'tO 

L ^ slaved disfcw* between s^-earTV, 

vx A i-hfanoTv^CQ ) otxrs 

The numerical approximation used in the program is 

AsO^t) it 

= expected incoming radiation. 

where QwOuOj is the angle between the reference area longitude and 
the sun projection on the equator. 

The quantity is incremented in 20 equal angular increments from sunset 
to sunrise. 
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Similarly 


A 10 

vo $ 0 ,<= p ,0 r r . s 0 \_ k ;) ao -^/^ 0 
j= i J 

The quantity r(t,) is an extrapolation (from NIMBUS II data) of the 
observed quantity throughout "all" t, angles. 
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Appendix F 


TABLES USED IN CALCULATION OF OUTGOING SOLAR RADIATION 


The two basic tables used in the albedo computations are presented with 


three independent variable breakpoints. Table F-l is used to determine 

r 1 


a value for the term 4- — — - L depending on the value of cos „ 

r If ( % ci S) i 

Table F-2 is used to determine a value for the term J -r : — r de* 

) j 

pending on the values of cos <Jj , sin P and cos ^ respectively. 


CoSC 

f (S~o) 

Ocean 

Land 

1.00 

1.0 

1.0 

0.954 

1.0 

1.0 

0.901 

1.0 

1.0 

0.848 

1.0 

1.0 

0.795 

1.0 

1.05 

0.742 

1.03 

1.1 

0.689 

1.1 

1.14 

0.636 

1.2 

1.18 

0.583 

1.3 

1.22 

0.530 | 

1.4 

1.28 

| 0.477 I 

1.6 

1.32 

0.424 ! 

1.8 

1.38 

! 0.371 j 

2.0 

1.42 

0.318 

2.2 

1.48 

0.265 

i 2.5 

1.52 | 

0.212 

1 2.8 

1.55 i 

0.159 

3.1 

1.58 ! 

0.106 

3.4 

1.6 1 

0.053 

3.7 

i.6 : 

0.0 

| 

4.0 

1.6 

i 


Table F-l 

K£)/r(£ =C 0 Tabular ions 


The independent variable isos J and it is divided into 20 intervals in 

' TabUs F_1 and F ‘ 2 are derl '' ed £ro “ data con- 
ained in the RADBAL computer program developed for Raschke, Vonder Haar 
Pasternak and Bandeen. * 
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Table F-2 entries are stored in the order (i, j, k) with the first index 
moving most rapidly. The breakpoints for i, j, and k are as follows: 


: 

i 

COS yj 

i 

— t.&< cos^ -,86603 

2 

- .86603 < cosy< -.5 

3 

-.5 < cosy^C 0. 

4 

0. < cos^< .5000 

5 

.5 < cos^< .86603 

! 6 

,86603 < cosj^c 1. 

i 

: 

j 

sin & 

1 

0 . < sin©< .17365 

2 

j .17365 < sine?< .34202 

3 

.34202 < sin^ .5000 

4 

.5000 < sin^x:' .64279 

5 

.64279 < sin£K 1 

k 

cos £ 

1 

.8195 < cos £ < 1 

2 

.5000 < cos ^ < .8195 

3 

.1736 ^ cos £ < .5000 
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Xtf (e'>4»/V') j 


OCEAN 

LAND 

1.2 

1.2 

1.3 

1.4 

1.4 

1.4 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

0.8 

1.1 

1.3 

1.4 

1.4 

1.4 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

0.5 

0.9 

1.25 

1.4 

1.4 

1.4 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

0.5 

0.9 

1.25 

1.4 

1.4 

1.4 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

0.7 

1.05 

1.25 

1.25 

1.2 

1.2 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.2 

1.2 

1.2 

1.25 

1.25 

1.3 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.05 

1.2 

1.2 

1.2 

1.25 

1.35 

0.95 

0.95 

0.95 

0.93 

0.93 

0.93 

0.85 

1.15 

1.2 

1.2 

1.2 

1.35 

0.9 

0.93 

0.93 

0.93 

0.9 

0.9 

0.55 

0.98 

1.15 

1.1 

1.1 

1.25 

0.83 

0.88 

0.93 

0.9 

0.87 

0.87 

0 . 35 

0.95 

1.15 

1.08 

1.05 

1.1 

0.7 

0.8 

0.85 

0.85 

0.83 

0.83 i 

! 0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

| 1.0 

1.0 

1.0 

1.02 

1.04 

1.04 i 

i 0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.98 

0.98 

1.0 

1.0 

1.04 

1.07 ! 

1 0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.93 

0.93 

0.97 

1.0 

1.04 

1.04 ! 

| 0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.8 

0.86 

0.93 

0.97 

0.97 

0.97 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

1 

0.74 

0.86 

0.92 

0.87 

0.84 
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Tabulations 


NOTE: The table values are shown reading from left to right 3 and are to 

be interpreted as shown in the following example. 
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1.2 I 
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i ! 

1.3 I 

: 4 

1 

i 1 

1.4 
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1 

i 

1.4 1 

6 

1 

i 

1.4 f 

1 

2 

! i 

0.8 

2 

! 2 

i l 

! l.i 

3 

2 ! 

i I 

1.3 

4 

1 2 ; 

l 
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1.4 

1 6 

i 

; • ; 

i 5 1 

i 

| 3 

o! 7 
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